Bloch oscillations in one-dimensional spinor gas 
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A force applied to a spin-flipped particle in a one-dimensional spinor gas may lead to Bloch 
oscillations of particle's position and velocity. The existence of Bloch oscillations crucially depends 
on the viscous friction force exerted by the rest of the gas on the spin excitation. We evaluate the 
friction in terms of the quantum fluid parameters. In particular, we show that the friction is absent 
for integrable cases, such as SU{2) symmetric gas of bosons or fermions. For small deviations from 
the exact integrability the friction is very weak, opening the possibility to observe Bloch oscillations. 
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Dynamics of ultracold atomic gases with internal 
(spinor) degree of freedom has been a focus of a number 
of recent experiments The observed collective phe- 
nomena have revived interest in earlier theoretical works 
Q on spin waves in helium and opened a possibility to 
study nonequilibrium dynamics of quantum liquids. 

Due to an unprecedented degree of experimental con- 
trol it is possible to excite a few atoms into a different 
hypcrfine internal state [1]. This leads to effective spin 
excitations which may be regarded as impurities mov- 
ing through the quantum liquid formed by the majority 
spins. A similar setup was investigated in the context of 
He^ and He^ mixtures 0, [1| . It was realized that an ex- 
ternal particle is " dressed" to form a collective excitation 
which energy-momentum relation at small momenta P is 
quadratic £{P) ^ iJ-d + P^/2M*. The correlations mani- 
fest themselves in quasiparticle effective mass M* being 
different from the bare mass M, as well as in the friction 
exerted on the quasiparticle by the rest of the liquid. 

The collective nature of the excitations is especially ap- 
parent in one-dimensional (ID) systems where the strong 
effects of interactions beyond mean field were recently ob- 
served in experiments with cold atoms Q . In addition to 
the strong mass renormalization Q , power law behavior 
of responce functions 0, 0, [l^l, the dispersion relation 
of the excitations £{P) was shown [ll| to be strongly 
modified by the interactions: parabolic at small P, it is 
actually a periodic function of the momentum with the 
period 2Tihn, see Fig.[T] Here n is a ID density of the gas 
and periodicity stems from the fact that total momentum 
Ptot — 2T:hn can be transferred to the g whole at 

no energy cost in thermodynamic limit. The periodicity 
of the dispersion relation drastically affects the dynam- 
ics of spin excitations under an influence of the external 
gravitational force F, which becomes uncompensated if 
the hypcrfine state of impurity atoms is insensitive to 
the magnetic field of the trap Indeed, momentum 

of impurity evolves according to P = F and its veloc- 
ity V = de{P)/dP is a periodic function of momentum 
and thus exhibits Bloch oscillations. Bloch oscillations 
are usually associated with an accelerated quantum par- 




FIG. 1: Dispersion relation of a spin excitation in ID Bose 
liquid. Inset: two-phonon processes leading to a dissipation. 
The arrows, with the slope given by the sound velocity c, rep- 
resent absorption and emission of long wavelength phonons. 



tide in presence of a static periodic potential, see, e.g., 
Ref. [ist for a recent experimental realization. However, 
the existence of Bloch oscillations in ID spinor conden- 
sates does not rely on a presence of an external periodic 
potential. It is the ID quantum liquid itself that pro- 
vides a quasi-periodic potential with the lattice spacing 
and thus 2Trhn reciprocal vector. 
Contrary to a static periodic potential, a ID quantum 
liquid exhibits quantum and thermal fluctuations. Be- 
cause of the latter the spin excitations with periodic dis- 
persion relation arc also subject to a dissipation. A pos- 
sibility to observe Bloch oscillations depends crucially on 
the strength of such a dissipation, i.e. the friction force 
exerted on the spin-flipped particle by the quantum liq- 
uid. Indeed, in a presence of friction the equations of 
motion for the spin excitation take the form 



P = F- kV: 



V = X ^ de/dP, 



(1) 



where k = k{V) is the coefficient of viscous friction 
[l^ . Consider, e.g., the strong coupling limit where 
e{P) + {2h^n^/M*) sm^[P/2hn) [ll and k w const. 
Integrating Eqs. ([T]), one finds for the drift velocity: 
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where the critical force is Fc = nhn/M*. As a result 
the drift velocity is a non-monotonic function of the ap- 
plied force and friction k. As discussed below, the lat- 
ter depends strongly on the temperature and interaction 
parameters, leading to a non-monotonic dependence of 
Vdrift on them too. In the strong forcing regime, F > Fc, 
the drift motion is superimposed with the Bloch oscil- 
lations with the period AT = 2nhn/y/F'^ - F^. Notice 
that both the period and the amplitude of the Bloch os- 
cillations diverge as F approaches Fc frome above. 

To find whether the Bloch oscillations regime is ex- 
perimentally accessible, one needs to have a detailed un- 
derstanding of the friction k{V). The friction originates 
from the interactions between the spin-excitation and the 
density excitations (phonons) of the rest of the gas. Since 
the latter propagate with the speed of sound c typically 
much larger than velocity of the particle V, the singlc- 
phonon processes do not lead to dissipation according 
to the Landau criterion [^. However, the two-phonon 
processes (see inset in Fig. [Ij do lead to energy and mo- 
mentum relaxation which scales as in 3D [1, 0] and 
as in ID [l^ (at small temperature and V < c). In 
this paper we evaluate the magnitude of friction and 
express it in terms of "Bose liquid" parameters. 

Dimensional arguments [l^ lead to k ~ T'^ / [c^fr'-n?), 
which for typical experimental parameters would make 
the critical force Fc rather large (see below). Fortu- 
nately the actual friction may be parametrically smaller 
than this estimate. The reason is that for certain sets of 
parameters the corresponding ID models are integrahle. 
One example is provided by SU (2) symmetric spinor liq- 
uid of bosons or fermions with point-like interactions [l6| , 
where the bare masses and interaction constants of both 
spin species are the same: M = m and G = g [l7| . 
The other example is the Tonks gas g ^ oo with M = m 
and an arbitrary interaction constant G between the spin 

As we explicitly ver- 



excitation and the liquid 
ify here in such integrablc cases the dissipation is absent 
[2 of and the spin excitations have infinite life-time even 
at finite temperature. Thus in the vicinity of integrable 
points the dissipation is strongly suppressed (see, e.g., 



Rcf. [2lj for experimental evidence in the spinless case). 



In a generic case we found that the friction coefhcient 
may be expressed through the thermodynamic properties 
of the gas. To this end one needs to know the chemical 
potentials of the majority /i(n) and minority fJ-d{n) spin 
species as functions of majority spins concentration n 
(the minority concentration is vanishingly small). For 
the friction coefficient at small velocity, ^ c, we found 



where the sound velocity is related to the compressibility 
as mc^ /n = dfi/dn and a ~ ii/dn^ , while 
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In the SU (2) symmetric case the equality of bare masses 
and interaction constants imply (5; = and ad = a. As 
a result the friction is absent k = 0. The same is true 
for the M ~ m Tonks gas, which is equivalent to an 
impurity moving through the non-interacting Fermi gas 
[231 . Finding the chemical potentials and the effective 
mass for non-integrable cases is a difficult task. Analytic 
progress may be achieved in two limits: 

Weak coupling limit: g <^ h?n/m and G fi^n/M. 
According to Bogoliubov theory the mean-field equations 
of state are ^ = gn and /i^ = Gn, while M* ~ M . Ex- 
tracting parameters as in Eq. ^ for both spin excitation 
and the liquid yields 
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Keeping the velocity dependence of the friction force, one 
finds k{V) = k(0)/(1 - V^/c^). As a result, the drift 
velocity of the spin excitations cannot surpass the speed 
of sound, Vdrift < c, making it hard to observe Bloch 
oscillations in weakly interacting gases. 

Strong coupling limit: (i) M ^ m and g,G oo. Both 
majority and minority spins may be considered as two 
free Fermi gases separated by impermeable mobile wall. 
Demanding the equal pressure of the two gases, one may 
find their total kinetic energy and thus respective chem- 
ical potentials: /i = {nhn)^ / {2m,) and Hd = {m/Aiy^^ii. 
Substituting it in Eqs. Q and employing that in this 
regime the effective mass is large M* /Ad ~ GM/h^n, one 
finds 
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(ii) M ~ m and g,G ^ h^n/m. In this case we treat the 
corresponding Fermi gases as being weakly interacting 
and find leading perturbative corrections to the chemi- 
cal potentials: /x = {■k'^H^ /2m){'n? — Wfi^n^ /3gm), while 
= {n^h'^/2m){n^ ~ Ari^n^il/g + l/3G]/m). Inserting 
it in the Eq. ([3]), one finds that in this order of expan- 
sion all the terms cancel each other. Such a cancelation 
is not expected if the order is kept in the chemical 
potentials. One can thus estimate the friction coefficient 
as 
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Since the effective mass of the spin excitations in the 
strong coupling regime is large M* ^ M Q , its velocity 
is small <C c and the friction coefficient is practically 



velocity-independent. It is in this regime where the Bloch 
oscillations are most likely to be observed. 

We turn now to the derivation of Eq. ([3]) . If the tem- 
perature is less than the chemical potential T ^ fi the 
density excitations of the majority spin gas may be de- 
scribed [i^l by the effective ID hydrodynamic Hamilto- 
nian 



77ph = I dx 
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hereafter H= 1. Here p{x) is the operator of density fluc- 
tuations on top of the uniform density n and its canoni- 
cally conjugate phase operator i9(a;) is related to the su- 
perfluid velocity Vs = Ox'^/m. The quadratic (Luttinger 
liquid) part of this Hamiltonian describes phonons with 
the linear dispersion relation uj{q) = c\q\. To have a con- 
sistent description of the interactions between phonons 
and the spin excitations, one needs to take into account 
non-linear interactions of phonons between themselves 
0. They are described by the terms ^ p{dx'&Y ^-nd 
~ p^ . The coefficients in front of them are dictated by 
Galilean invariance for the former, and by expansion of 
p{n + p{x)) up to the second order in p for the latter. 

The spin excitation may be thought of as a quantum 
particle described by the canonically conjugated coor- 
dinate X and momentum P. Its interactions with the 
density fluctuations are encoded in the chemical poten- 
tial pd{n + p{X)) « (mcVn)(l - 5i)p{X) + {ad/2)p'{X), 
cf. Eq. ([3]). The interactions of this particle with the 
superfluid velocity Vs{X) may be found Q by notic- 
ing that in the reference frame where Vs = ^ the en- 
ergy of the particle with momentum P is given by the 
dispersion relation s{P). In the laboratory frame mo- 
mentum of such particle is P + Mvs, while its energy 
is Ey^iP + Mvs) = e{P) + Pi's + Mv1/2. Changing 
to a momentum in the laboratory frame, and keeping 
terms up to the second power in Vg^ one finds ey^{P) = 
£{P) + {P - MV)vs + {APe"{P) - M)vl/2, where the 
particle velocity \&V = s'{P). At a sufficiently small ve- 
locity V = P/M* and £"(0) = 1/M* one thus finds for 
the Hamiltonian of the spin excitations 



{P + 5Mvsf 
2M* 



5Mi 



-(1 



^/)p+^P^(9) 



where 5M = M* — M . In the absence of interactions 
5M = M* ~ M = and Hamiltonian ^ is independent 
of the fluid velocity. The terms containing p = p{X) and 
Vs — dx'd{X)/m can be regarded as effective interaction 
potential dependent on the particle coordinate X. This 
introduces a preferential frame for the moving particle. 

It is convenient to perform canonical transformation 
of the particle momentum P -\- SMvg — > P along with 
the fluid "coordinate" p{x) — [SM / m) pd{x) — > p{x), 
where pd{x) = 5{x — X) is the density of the particle. 
The changes induced by this transformation to the fluid 






FIG. 2: Second order diagrams contributing to 2-phonon am- 
plitudes. Spin excitation is represented by a full line, while 
phonons by wavy lines. Diagrams a,b and c contribute to Fp, 
Eq. p2|l . while diagram d represents F^, Eq. (|13p . 



Hamiltonian, Eq. ([H]), are absorbed in the modified im- 
purity Hamiltonian 
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(10) 

The second term here describes processes in which one 
phonon is absorbed or emitted, while the third one is 
responsible for the two-phonon processes. Due to the 
quadratic dispersion relation of the spin excitation and 
linear dispersion of phonons, the one-phonon processes do 
not lead to dissipation, which is just another statement 
of Landau criterion. 

We thus focus on the two-phonon amplitude. The lat- 
ter originates from the last term in Eq. (fTO|) as well as 
the following second order processes: (i) second order 
in p{X) — Jdxp{x)pd{x) interaction vertex. Fig. [2^,b; 
(ii) first order in ppd and first order in p'^, Fig. [2t or 
p{dx^)^ phonon nonlinearity vertexes, Fig. [^Ji. It is the 
destructive interference of these second order processes 
which is responsible for the partial or even complete sup- 
pression of the dissipation. Evaluating the corresponding 
diagrams according to the standard rules, one derives an 
effective Hamiltonian of the spin excitation 24 1 



^~ljp[pix)T~lMdxm)]\ (11) 



where the effective two-phonon amplitudes are given by 
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The momentum relaxation rate may be evaluated in 
the second order in the two-phonon amplitudes Fp^^ , av- 
eraged over the Luttinger (quadratic) part of the phonon 
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Hamiltonian ([8]). This does not assume smallness of 
the ampUtudcs Fp ,?, but rather gives the leading low- 
temperature, T <^ /i, contribution. The semiclassical 
equation of motion for the spin excitation acquires a form 



P 



^gn(g.gv). (14) 



where n(q, lo) is the imaginary part of the Fourier trans- 
form of 9{t){[p'^{x,t), p'^{0,0)]) response function of the 
phonon gas. The latter is given by 



n(<z,^) 



-TT I coth 
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coth 
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(15) 

For spin excitations of small velocity, V <^ one finds 



n(g,gy) 



V 



T sinh'(cg/4r) 



(16) 



Substituting Eqs. ^ and ^ in Eq. (dH), one 

finds the friction force P = —kV with the friction co- 
efficient K given by Eq. Notice that in the inte- 
grable SU{2) symmetric case the two-phonon amplitudes 
Tp = M*c^ I'n? and F^ = M* jw? are finite, while the 
dissipation rate vanishes due to the interference of the 
density and current excitations, cf. Eq. (|14p . 

Finally we give some estimates for the critical force F^.. 
For ®^Rb gas with the linear density n = lO^cm"^ and 
the ID interaction constant g of the same order as h^n/m, 
one estimates the sound velocity as c « \.Ahn/m ~ 
Icm/s. Taking M = m and M* w 1.3 m 0], and 
the temperature T « 0.5 mc^ « 5 x 10^' K, one finds 
for the naive critical force Fi"' « n^T* / {c^fi^nM*) « 
4 X 10~^^ N. For comparison, the gravitational force act- 
ing on s'^Rb atom is 1.4 x lO^^'^N, i.e. about 300 times 
weaker. However, closeness to the integrability with say 
G/g = 1.03 decreases the actoa/ critical force Fc by about 
three orders of magnitude, see Eqs. ([5]) or ([7]) . This makes 
the critical force to be well below the gravitational one, 
making it possible to observe Bloch oscillations of the 
spin-flipped particle falling in the gravitational field. 
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